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Abstract 

The production of ion beams from the interaction of a circularly polarized laser pulse with a nanometric 
double-layer target is discussed in the regime where all electrons are expelled from the target by laser radi- 
ation pressure. Quasi-monochromatic, well-collimated ion beams are observed in two-dimensional particle- 
in-cell simulations. The ion beam properties are derived from a simple analytical model, and the possibility 
to control those properties by using a laser-pulse with sharp-rising edge is discussed. Potential application 
to hadron-therapy is finally considered. 



1. Introduction 

Generating high-energy ion beams by irradiating 
a solid target with an ultra-intense laser pulse has 
many potential applications from medicine to con- 
trolled nuclear fusion While different mecha- 
nisms of ion acceleration have been demonstrated, 
ion acceleration at currently available laser inten- 
sities mainly follows from the plasma expansion 
driven by the hot, laser-created, electrons [2J. Ex- 
perimentally, protons with energy up to 70 MeV 
have already been produced [3J. However, con- 
trolling their energy distribution remains a crucial 
problem and alternative acceleration mechanisms 
are considered to overcome this issue 0, 0, d, 0| • 

In this paper, we revisit ion acceleration from a 
double-layer target and propose a new laser-based 
acceleration scheme that allows an excellent con- 
trol of the ion beam properties. This mechanism 
relies on the complete expulsion of electrons from 
an ultra-thin double-layer target by laser radia- 
tion pressure. The target is almost instantaneously 
transformed into a capacitor-like structure, where 
the ions of the first layer and the electron cloud 
play the role of the anode and cathode, respec- 
tively. A strong, quasi-homogeneous, electrostatic 
field (~ lOOTV/m) is built up between the two 



* Corresponding author 
Email address: mickael.grechagmail.com (M. Grcch) 
Preprint submitted to Elsevier 



layers. The ions of the second layer, with a larger 
charge-over-mass ratio, are then accelerated in a 
way similar to ions in a conventional linear accel- 
erator. Later on, once electrons escape from the 
ion layers, light ions can gain further energy in the 
electrostatic field of the bare first ion layer. 

Two-dimensional (2D) particle-in-cell (PIC) sim- 
ulations and an analytical model show that high- 
energy (> 100 MeV) ion beams can be generated 
over a few laser periods only. These beams are char- 
acterized by small and adjustable energy (< 10 %) 
and angular (< 5°) dispersions, which makes them 
suitable for applications such as hadron-therapy. 

2. Ion beams from a double-layer target 

2.1. Basic picture of ion acceleration 

Double-layer targets for laser-based ion acceler- 
ation have recently attracted a lot of attention as 
they allow to control the energy dispersion of the 
accelerated ions 0, The first target layer is 
usually made of heavy ions with mass mh ^> ni p 
(m p = 1836 is the proton mass) and charge Zh. 
The second target layer (or a small dot placed at 
the rear-side of the first layer) consists of ions with 
a larger charge-over-mass ratio (Zi/mi > Zh/rrih). 
In what follows, we denote by nh, dh and Wh the 
atomic density, thickness and transverse size of the 
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first layer, respectively. Similarly m, di and wi de- 
note the corresponding parameters for the second 
target layer or structured dot (if wi < Wh)- For 
the sake of simplicity all quantities in this paper 
are expressed in normalized units: masses, charges, 
densities and distances are normalized to the elec- 
tron mass m e , electron charge e, critical density 
n c = eo m e oj L /e 2 (eo is the permittivity of vacuum 
and lul is the laser frequency) and inverse laser wave 
number k^ 1 = c/lul (c is the light velocity), respec- 
tively In addition, the laser field amplitude a^,, as 
well as all other electric fields, are given in units of 
the Compton field Ec = m e ojL c/e. 

When irradiated by an intense laser pulse (at 3> 
1), the target (or a part of the target) is quickly ion- 
ized. Electrons are expelled from the target, and a 
strong electrostatic field is built-up. Light ions of 
the second layer are accelerated in this field. The 
condition for generation of quasi-monochromatic, 
well-collimated, light ion beams is simply that 
the accelerating electrostatic field remains quasi- 
homogeneous over the second target layer during 
the acceleration process [H, I n this paper, we 
revisit ion acceleration from a double-layer target 
in the regime where the laser radiation pressure ef- 
fectively expels all electrons from the target. 

2.2. 2D PIC simulations 

Two-dimensional simulations are performed us- 
ing the PIC code C ALDER 0. A laser pulse with 
field amplitude = 100 is focused at normal in- 
cidence on an ultra-thin (dh,di <C 1) double- layer 
target. A circularly polarized laser pulse with flat- 
top temporal intensity profile is used in order to 
make electron sweepin g by the laser radiation pres- 
sure more efficient llOj. The pulse duration is 
t p ~ 62.8 (10 optical cycles) only. Along the trans- 
verse y— direction, the laser intensity follows a 6-th 
order super-Gaussian profile. Two values for the 
laser focal-spot full-width at half-maximum have 
been considered: wl = 62.8 and 125.6 (10 and 20 
laser wavelengths, respectively). The first target 
layer is made of fully ionized carbon (Z^ = 6, rrih = 
12 m p ) with atomic density = 58 and thickness 
dh = 0.08. Its transverse size is Wh ~ 200 (~ 30 
laser wavelengths) , which is also the transverse size 
of the simulation domain. A likewise fully ionized 
hydrogen dot (Zi = 1, mj = m p ) with transverse 
size wi = 6.28 (one laser wavelength) and thickness 
di = 0.16 is placed at the rear-side of the first layer. 
Two atomic densities, ni = 6 and ni = 12, have 



been considered for the hydrogen dot. The simula- 
tions results are summarized in Figs. [1] and [5] 

Four laser-periods after the beginning of the in- 
teraction, all electrons in a region with diameter 
tui ~ 70 ~ wl irradiated by the laser pulse have 
been removed from the target (Fig. [TJi). In con- 
trast to previous works 0, Q where electrons are 
strongly heated and generate an accelerating sheath 
at the target rear-side, electrons here form a com- 
pact bunch that propagates with the front of the 
laser pulse. The target then exhibits a capacitor- 
like structure and a quasi-homogeneous electro- 
static field is built up between the bare heavy ion 
layer and the electron bunch (Fig. QJ1) . Light ions 
are accelerated in this electrostatic field and can 
reach high energies in a few optical cycles only. In 
Fig. [2h>, the light ion energy e is shown as a func- 
tion of time. Four optical cycles after the beginning 
of the interaction, protons have already reached 
e ~ 70MeV (blue line, square). 

Due to the quasi-homogencity of the accelerating 
field, both relative energy dispersion Ae/e ~ 10% 
and aperture angle Ad < 5° remain small (Figs, 
and[2t). Moreover, both of them are, at this stage 
of the acceleration process, sensitive to the charge 
density in the proton dot: the electrostatic repul- 
sion between the protons themselves is therefore the 
main source of energy and angular dispersion. 

This early stage of proton acceleration is main- 
tained as long as the distance between the electron 
cloud and the bare heavy ion layer remains small 
compared to their transverse size w±. At the end 
of this stage (at t ~ 50 for the smaller laser focal 
spot and t ~ 80 for the larger one), protons with 
energies of the order of 200 MeV and 480 MeV, re- 
spectively, have been produced. Later on, the elec- 
trostatic field seen by the protons strongly decreases 
due to geometrical effects (Figs.[Tb,c,e,f and[5h.) and 
the acceleration process slows down (Fig. Wp)- 

Figure [3Jd also shows that, while the proton en- 
ergy in the early stage of ion acceleration does not 
depend on the transverse width wl of the laser focal 
spot, increasing wl allows to enhance the final pro- 
ton energy by delaying multi-dimensional effects. 
Furthermore, this later stage is associated with a 
slight increase of both the energy and angular dis- 
persion of the proton bunch (Figs.[2fc and[2jl). This 
is the sign that both dispersions result from elec- 
trostatic repulsion in the proton dot itself, which 
evolves on a time scale different from that of the 
acceleration process. Nevertheless, for the lower hy- 
drogen density (n; = 6), energy and angular disper- 
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Figure 1: (a,b,c) Snapshots of the electron density (gray scale), carbon contour plot (blue) and proton contour plot (red), 
and (d,e,f) corresponding electrostatic field distribution at three different times after the beginning of the interaction: (a,d) 
4 optical cycles; (b,e) 14 optical cycles and (c,f) 30 optical cycles. The colormap refers to panels (d-f). In this simulation, 
wl = 62.8 (10 laser wavelengths) and n; = 6. 



sion seem to saturate at the end of the simulation to 
values not exceeding 10 % and 3° , respectively. Fi- 
nally, we point out that while Ae/e is only slightly 
sensitive to the transverse size of the laser pulse, A9 
is decreased for the larger focal spot. This prompts 
us to suggest that the transverse inhomogeneity of 
the accelerating field also plays a role in the proton 
angular dispersion. 

3. Analytical model 

On the contrary to previous studies [1, the 
acceleration mechanism discussed here relies on the 
complete and directed expulsion of all electrons 
from the double-layer target. This occurs if the 
laser radiation pressure ~ a\ is sufficient to over- 
come both the electron thermal pressure, and the 
maximal electrostatic pressure a 2 ot /2 resulting from 
complete electron-ion separation. The electrostatic 
field a to t = a h + a t , where a h = Z h n h d h and 
ai = Zi m di , is the electrostatic field between the 
bare ion layers and the electron cloud. Using a 
circularly polarized laser pulse allows to maintain 
the electron temperature to a rather low level 
so that one obtains the simple condition for com- 
plete electron removal a\ > a 2 ot /2. This condi- 
tion defines a threshold for the laser field amplitude 



o-l > a>tot beyond which electrons are removed from 
the target and propagate as a compact bunch at the 
front of the laser pulse [ToT ]. 

3.1. Linear plasma acceleration 

Previous numerical simulations show that early 
light ion acceleration (typically for t < wl) is one- 
dimensional (ID). It occurs in a monotonically in- 
creasing accelerating field that preserves the order 
of the different ion sheets that make up the light ion 
layer. Light ions therefore see a constant acceler- 
ating field whose amplitude depends only on their 
initial position. The light ion mean energy e and 
energy dispersion Ae (in units of m e c 2 ) can thus 



be computed exactly [111 ]: 



Ae/i„PA 



mi 
mi 



(1) 



^l+(l + ai /a h ) 2 t 2 /t 2 



(2) 



where the index linPA accounts for linear plasma 
acceleration (with obvious analogy to conventional 
acceleration techniques), and t r = mi/(Zi<Xh) de- 
notes the time over which light ions gain relativis- 
tic energies. For f<t r , they behave classically and 
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Figure 2: Temporal evolution of: (a) the accelerating elec- 
trostatic field a e . a . at the position of the proton bunch; 
(b) the proton energy e; (c) the energy dispersion Ae/e; 
(d) the angular aperture A9. For n; = 6, wl = 10 Xl (blue, 
square), n; = 12, = 10 (green, triangle) and n; = 6, 
wl = 20 Xl (black, circle). Dashed curves account for ana- 
lytical predictions of Sec. [3] Note that in panels (a) and (b) 
the blue and green curves are superimposed. 



the mean energy and relative energy dispersion are 
e<» ~ m l t 2 /{2t 2 r ) and (Ae/e)( c > - 2ai/a h , respec- 
tively. If light ions experience the field ~ ah over 
a time t 3> t r , they gain ultra-relativistic energies 



mi t/t r , and their relative energy dispersion 



(•ur) 

simply reads (Ae/e)( ur ) ~ ai/au- 

These analytical estimates are compared to the 
2D simulation results of Sec. 12.21 In the first stage 
of the acceleration process, light ions see a constant 
accelerating field a to t ~ 28 (Fig. [2^,). The result- 
ing temporal evolution of the ion energy e as pre- 
dicted by Eq. (fTJ) is in excellent agreement with the 
numerical results during the first acceleration stage 
(Figs.O)). A good agreement is also found between 
the predicted value of Ae/e and numerical results 
(Fig. . In particular, Ae/e cx ai/ah is kept small 
by ensuring that the accelerating field ah remains 
large compared to a; , which determines the ratio of 
the areal charges of both ion layers, and limits the 
brightness of the resulting ion source. 



3.2. Long time behavior and DCE 

The ID behavior of light ion acceleration breaks 
down once the electron cloud is separated from the 
heavy ion layers by a distance of the order of its 
transverse size w± ~ w L . Because electrons leave 
the target at a velocity close to c, this defines a 
characteristic time timPA ~ Wjl ~ wl over which 
the capacitor-like acceleration is maintained. On 
longer times, electrons eventually escape from the 
ion layers, and their contribution to the later stage 
of the ion acceleration process can be neglected. 
Light ions thus can gain further energy in the elec- 
trostatic field due to the bare heavy ion layer in a 
way similar to the so-called directed Coulomb ex- 
plosion (DCE) Q. Assuming Zh/nih <C 1, the en- 
ergy gained during this later stage can be derived 
from the conservation of the total light ion energy: 



£dcb ~ Zi a h w±/4— eunPA/2 . 



(3) 



The first term in the right-hand-side of Eq. ([3]) is 
similar to estimates obtained in Ref. [4[ . The second 
term accounts for the fact that, in contrast to what 
occurs in "classical" DCE, part of the light ions 
potential energy has already been transformed into 
kinetic energy in the first acceleration stage. 

As suggested by our numerical simulations, en- 
ergy dispersion during this later stage mainly fol- 
lows from the electrostatic repulsion between the 
accelerated light ions. A first attempt to estimate 
this effect is given in Ref. [ll|, however, the method 
used there is not correct. Here, we propose to esti- 
mate the energy dispersion by considering the ion 
bunch at the end of its acceleration stage as an ho- 
mogeneously charged sphere with radius R ~ wi/2, 
expanding due to its own charge Qi (in units of 
en c /k L ). The final energy dispersion can then be 
derived as a function of the final ion energy e: 



Ae_oc-E 



(4) 



where e p = ZiQi/ (2nwi) is the initial potential en- 
ergy of a light ion on the outer shell of the sphere. 

3.3. Scaling laws 

A question arises from previous calculations: 
which of the two accelerating stages provides the 
main contribution to the light ion acceleration? By 
comparing Eqs. (Q]) and ©, one obtains that ions 
with the final energy e < mi /8 gain most of their en- 
ergy during the later, DCE-like, acceleration stage. 
For these ions, the final energy is e ~ ZiahW^/A 
and energy dispersion Ae/e ~ 4 y^e p /e follows 



mainly from Coulomb repulsion. Both quantities 
can be controlled by adjusting the target parame- 
ters. Moreover, since w± ~ Wl, the condition for 
electron removal, ah > a>hi defines a threshold for 
the laser power and we obtain that the light ion 
energy scales as e oc Z\ t/Pl/4. 

In contrast to that, more energetic ions (in par- 
ticular relativistic ions) quickly separate from the 
heavy ion layer so that they gain most of their en- 
ergy during the early, linPA, stage. Using tu n pA ~ 
w± ~ Wi and cll > dhi Eq. JT]) can be rewritten 
as a function of the laser power. One then obtains 
a threshold power, Pl > mf/Zf, above which rela- 
tivistic ions with energy e oc Z\ V ' Pl are expected. 

4. Application to hadron-therapy 

Hadron-therapy is certainly one of the most ap- 
pealing applications of laser-created ion beams. It 
typically requires ions with energy of the order of 
200 McV/nucleon. For such energies, mi/8 < e < 
mi, both stages of the acceleration process are im- 
portant (see Sec. I3.3|) . The final energy and energy 
dispersion arc thus obtained from Eqs. |T|) and Q 
and Eqs. (0) and ([4]), respectively We now give 
the characteristic parameters of a target designed 
for generation of high-quality proton or carbon ion 
beams. For given ion energy (let us say 200 MeV 
per nucleon) and charge-over-mass ratio, the quan- 
tity ah w± is fixed. Considering w± ~ 62.8 (10 laser 
wavelengths), we obtain that accelerating field am- 
plitudes ah ~ 16 and ah ~ 32 are required to cre- 
ate the proton and carbon ion beams, respectively 
For proton acceleration, a nanomctric (~ 5 nm- 
thick considering a 1 jum-lascr wavelength), fully- 
ionized, diamond-like carbon foil at solid density 
can be used as a first layer to generate such an ac- 
celerating field. For carbon ion beam generation, 
a gold target is here considered. Assuming partial 
ionization of the gold atoms (Zh ~ 50), the required 
target thickness is ~ 1 nm. The second ion layer 
transverse width wi ~ 12.6 (two laser wavelengths) 
is chosen <C w± in order to minimize the transverse 
inhomogeneity of the accelerating field. Consider- 
ing Ae/e ~ 2%, we finally obtain from Eqs. (J2J) 
and (|H) that Coulomb repulsion between the ac- 
celerated ions limits the accelerated charge Qi to 
~ 0.4 pC for protons and ~ 0.9 pC for carbon ions 
per shot. We would like to stress that this limita- 
tion is generic to short, non- neutral ion bunches. 

In order to ensure that the laser radiation pres- 
sure is sufficient to remove all electrons from the 



target, laser intensities above 3.5 x 10 20 W/cm 2 and 
1.4 x 10 21 W/cm 2 (correspondingly laser powers ex- 
ceeding 360 TW and 1.5 PW, and laser energies of 
a few tens of J) are required depending on whether 
one considers proton or carbon ion acceleration. To 
accelerate a sufficient number of ions, repetition 
rates in the kHz-domain are necessary. 

5. Conclusion 

In summary, the generation of quasi- 
monochromatic, well-collimated, light ion beams 
from a nanometric double-layer target is discussed 
in the regime where the laser radiation pressure is 
sufficient to remove all electrons from the target. 
The resulting ion acceleration, which is shown 
to occur in two steps, is described using 2D PIC 
simulations and an analytical model. Both, the 
light ion source energy and energy dispersion are 
shown to depend mainly on the target design. On 
the contrary to previous findings, we show that the 
light ion source energy dispersion depends mainly 
on the electrostatic repulsion in the accelerated 
bunch itself, which limits the source brightness. 
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